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Abstract We review a recent proposal of a first principles approach to the electronic struc- 
ture of materials with strong electronic correlations. The scheme combines the 
GW method with dynamical mean field theory, which enables one to treat strong 
interaction effects. It allows for a parameter-free description of Coulomb in- 
teractions and screening, and thus avoids the conceptual problems inherent to 
conventional "LDA+DMFT", such as Hubbard interaction parameters and dou- 
ble counting terms. We describe the application of a simplified version of the 
approach to the electronic structure of nickel yielding encouraging results. Fi- 
nally, open questions and further perspectives for the development of the scheme 
are discussed. 

Keywords: Strongly correlated electron materials, first principles description of magnetism, 
dynamical mean field theory, GW approximation 



Introduction 

The development of density functional theory in combination with rapidly 
increasing computing power has led to revolutionary progress in electronic 
structure theory over the last 40 years. Nowadays, calculations for materials 
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with large unit cells are feasible, and applications to biological systems and 
complex materials of high technological importance are within reach. 

The situation is less favorable, however, for so-called strongly correlated 
materials, where strong localization and Coulomb interaction effects cause 
density functional theory within the local density approximation (DFT-LDA) 
to fail. These are typically materials with partially filled d- or f-shells. Fail- 
ures of DFT-LDA reach from missing satellite structures in the spectra (e.g. in 
transition metals) over qualitatively wrong descriptions of spectral properties 
(e.g. certain transition metal oxides) to severe qualitative errors in the calcu- 
lated equilibrium lattice structures (e.g. the absence of the Ce a — 7 transition 
in LDA or the 30% error on the volume of <5-Pu). While the former two situa- 
tions may be at least partially blamed on the use of a ground state theory in the 
forbidden range of excited states properties, in the latter cases one faces a clear 
deficiency of the LDA. We are thus in the puzzling situation of being able to 
describe certain very complex materials, heading for a first principles descrip- 
tion of biological systems, while not having successfully dealt with others that 
have much simpler structures but resist an LDA treatment due to a particular 
challenging electronic structure. The fact that many of these strongly corre- 
lated materials present unusual magnetic, optical or transport properties has 
given additional motivation to design electronic structure methods beyond the 
LDA. 

Both, LDA+U [1-4] and LDA+DMFT [5-8] techniques are based on an 
Hamiltonian that explicitly corrects the LDA description by corrections stem- 
ming from local Coulomb interactions of partially localized electrons. This 
Hamiltonian is then solved within a static or a dynamical mean field approxi- 
mation in LDA+U or LDA+DMFT respectively. In a number of magnetically 
or orbitally ordered insulators the LDA underestimation of the gap is success- 
fully corrected by LDA+U (e.g. late transition metal oxides and some rare earth 
compounds); however its description of low energy properties is too crude to 
describe correlated metals, where the dynamical character of the mean field is 
crucial and LDA+DMFT thus more successful. Common to both approaches 
is the need to determine the Coulomb parameters from independent (e.g. "con- 
strained LDA") calculations or to fit them to experiments. Neither of them thus 
describes long-range Coulomb interactions and the resulting screening from 
first principles. This has led to a recent proposal [9] of a first principles elec- 
tronic structure method, dubbed "GW+DMFT" that we review in this article. 
Similar advances have been presented in a model context [10]. 

The paper is organized as follows: in section 1 we give a short overview 
over the parent theories (GW, DMFT, and LDA+DMFT), while in section 2 we 
introduce a formal way of constructing approximations by means of a free en- 
ergy functional. The form of this functional defining the GW+DMFT scheme 
is discussed in section 3; section 4 presents different conceptual as well as tech- 
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nical issues related to this scheme. Finally we present results of a preliminary 
static implementation combining GW and DMFT, and conclude the paper with 
some remarks on further perspectives for the development of the GW+DMFT 
scheme. 

1. The parent theories 
The GW Approximation 

Even if density functional theory is strictly only applicable to ground state 
properties, band dispersions of sp-electron semi-conductors and insulators have 
been found to be surprisingly reliable - apart from a systematic underestima- 
tion of band gaps (e.g. by ~ 30% in Si and Ge). 

This underestimation of bandgaps has prompted a number of attempts at 
improving the LDA. Notable among these is the GW approximation (GWA), 
developed systematically by Hedin in the early sixties [11]. He showed that 
the self-energy can be formally expanded in powers of the screened interac- 
tion W, the lowest term being iGW, where G is the Green function. Due to 
computational difficulties, for a long time the applications of the GWA were 
restricted to the electron gas. With the rapid progress in computer power, ap- 
plications to realistic materials eventually became possible about two decades 
ago. Numerous applications to semiconductors and insulators reveal that in 
most cases the GWA [12-14] removes a large fraction of the LDA band-gap 
error. Applications to alkalis show band narrowing from the LDA values and 
account for more than half of the LDA error (although controversy about this 
issue still remains [15]). 

The GW approximation relies on Hedin's equations [11], which state for the 
self-energy that 

E(l,2) = -t J d3d4v(l,4)G(l,3) 6G 5 ^ 2) (1) 

where v is the bare Coulomb interaction, G is the Green function and <\> is an 
external time-dependent probing field. We have used the short-hand notation 
1 = [x\t\). From the equation of motion of the Green function 

G- 1 = i— - H - £ (2) 
at 

H = h + (f> + V H (3) 
ho is the kinetic energy and Vh is the Hartree potential. We then obtain 



«ty(4) -- d v- 2 )[ d v-v + -sm- 



(52(3,2) 
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= _ 5(3 _2) e -i(3,4)-fgf (4) 

where e 1 is the inverse dielectric matrix. The GWA is obtained by neglecting 
the vertex correction 5T,/54>, which is the last term in (4). This is just the 
random-phase approximation (RPA) for e _1 . This leads to 

E(l,2) = iG(l,2)W(l,2) (5) 
where we have defined the screened Coulomb interaction W by 



W(l,2) = J d3v(l,3)e _:l (3,2) 
The RPA dielectric function is given by 



(6) 



e = l-vP (7) 

where 



P(r,r';u) = -2i f — G(r, r'; uj + u')G(r', r; J) 
J 2tt 

OCC UTIOCC 

= 2]T E ^(r)^(r')^(r)^(r') (8) 

i 3 
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with the Green function constructed from a one-particle band structure {?/>«, Ei}. 
The factor of 2 arises from the sum over spin variables. In frequency space, 
the self-energy in the GWA takes the form 

S(r, r'; w) = ^- [ duj'e^' G(r, r'; oj + uj')W(r, r'; J) (10) 

2-7T J 

We have so far described the zero temperature formalism. For finite tempera- 
ture we have 

2 

P(r,r';«/ n ) = - ^ G( r , r'; iv n + iu k )G(r', r; iu k ) (11) 

Y>{r,r';iuj n ) = ~^G{r,r' ';ioj n + ii/ k )W(r,r' ';iv k ) (12) 

In the Green function language, the Fock exchange operator in the Hartree- 
Fock approximation (HFA) can be written as iGv. We may therefore regard 
the GWA as a generalization of the HFA, where the bare Coulomb interaction 
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v is replaced by a screened interaction W. We may also think of the GWA as a 
mapping to a polaron problem where the electrons are coupled to some bosonic 
excitations (e.g., plasmons) and the parameters in this model are obtained from 
first-principles calculations. 

The replacement of v by W is an important step in solids where screen- 
ing effects are generally rather large relative to exchange, especially in metals. 
For example, in the electron gas, within the GWA exchange and correlation 
are approximately equal in magnitude, to a large extent cancelling each other, 
modifying the free-electron dispersion slightly. But also in molecules, accu- 
rate calculations of the excitation spectrum cannot neglect the effects of corre- 
lations or screening. The GWA is physically sound because it is qualitatively 
correct in some limiting cases [19]. 

The success of the GWA in sp materials has prompted further applications 
to more strongly correlated systems. For this type of materials the GWA has 
been found to be less successful. Application to ferromagnetic nickel [16] il- 
lustrates some of the difficulties with the GWA. Starting from the LDA band 
structure, a one-iteration GW calculation does reproduce the photoemission 
quasiparticle band structure rather well, as compared with the LDA one where 
the 3d band width is too large by about 1 eV. However, the too large LDA ex- 
change splitting (0.6 eV compared with experimental value of 0.3 eV) remains 
essentially unchanged. Moreover, the famous 6 eV satellite, which is of course 
missing in the LDA, is not reproduced. These problems point to deficiencies in 
the GWA in describing short-range correlations since we expect that both ex- 
change splitting and satellite structure are influenced by on-site interactions. In 
the case of exchange splitting, long-range screening also plays a role in reduc- 
ing the HF value and the problem with the exchange splitting indicates a lack 
of spin-dependent interaction in the GWA: In the GWA the spin dependence 
only enters in G but not in W. 

The GWA rather successfully improves on the LDA errors on the bandgaps 
of Si, Ge, GaAs, ZnSe or Diamond, but for some materials, such as MgO and 
InN, a significant error still remains. The reason for the discrepancy has not 
been understood well. One possible explanation is that the result of the one- 
iteration GW calculation may depend on the starting one-particle band struc- 
ture, since the starting Green's function is usually constructed from the LDA 
Kohn-Sham orbitals and energies. For example, in the case of InN, the starting 
LDA band structure has no gap. This may produce a metal-like (over)screened 
interaction W which fails to open up a gap or yields a too small gap in the GW 
calculation. A similar behaviour is also found in the more extreme case of NiO, 
where a one-iteration GW calculation only yields a gap of about 1 eV starting 
from an LDA gap of 0.2 eV (the experimental gap is 4 eV) [17, 12]. This prob- 
lem may be circumvented by performing a partial self-consistent calculation in 
which the self-energy from the previous iteration at a given energy, such as the 
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Fermi energy of the centre of the band of interest, is used to construct a new set 
of one-particle orbitals. This procedure is continued to self-consistency such 
that the starting one-particle band structure gives zero self-energy correction 
[17, 12, 18]. In the case of NiO this procedure improves the band gap consid- 
erably to a self-consistent value of 5.5 eV and at the same time increases the 
LDA magnetic moment from 0.9 \i& to about 1.6 \ib much closer to the exper- 
imental value of 1.8 \xb ■ A more serious problem, however, is describing the 
charge-transfer character of the top of the valence band. Charge-transfer insu- 
lators are characterized by the presence of occupied and unoccupied 3d bands 
with the oxygen 2p band in between. The gap is then formed by the oxygen 2p 
and unoccupied 3d bands, unlike the gap in LDA, which is formed by the 3d 
states. A more appropriate interpretation is to say that the highest valence state 
is a charge-transfer state: During photoemission a hole is created in the tran- 
sition metal site but due to the strong 3d Coulomb repulsion it is energetically 
more favourable for the hole to hop to the oxygen site despite the cost in en- 
ergy transfer. A number of experimental data, notably 2p core photoemission 
resonance, suggest that the charge-transfer picture is more appropriate to de- 
scribe the electronic structure of transition metal oxides. The GWA, however, 
essentially still maintains the 3d character of the top of the valence band, as in 
the LDA, and misses the charge-transfer character dominated by the 2p oxygen 
hole. A more recent calculation using a more refined procedure of partial self- 
consistency has also confirmed these results [18]. The problem with the GWA 
appears to arise from inadequate account of short-range correlations, probably 
not properly treated in the random-phase approximation (RPA). As in nickel, 
the problem with the satellite arises again in NiO. Depending on the starting 
band structure, a satellite may be reproduced albeit at a too high energy. Thus 
there is a strong need for improving the short-range correlations in the GWA 
which may be achieved by using a suitable approach based on the dynamical 
mean-field theory described in the next section. 

Dynamical Mean Field Theory 

Dynamical mean field theory (DMFT) [20] has originally been developed 
within the context of models for correlated fermions on a lattice where it has 
proven very successful for determining the phase diagrams or for calculations 
of excited states properties. It is a non-perturbative method and as such ap- 
propriate for systems with any strength of the interaction. In recent years, 
combinations of DMFT with band structure theory, in particular Density func- 
tional theory with the local density approximation (LDA) have emerged [5, 6]. 
The idea is to correct for shortcomings of DFT-LDA due to strong Coulomb 
interactions and localization (or partial localization) phenomena that cause ef- 
fects very different from a homogeneous itinerant behaviour. Such signatures 
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of correlations are well-known in transition metal oxides or f-electron systems 
but are also present in several elemental transition metals. 

Combinations of DFT-LDA and DMFT, so-called "LDA+DMFT" techniques 
have so far been applied - with remarkable success - to transition metals (Fe, 
Ni, Mn) and their oxides (e.g. La/YTi0 3 , V 2 3 , Sr/CaV0 3 , Sr 2 Ru0 4 ) as well 
as elemental f-electron materials (Pu, Ce) and their compounds [7]. In the most 
general formulation, one starts from a many-body Hamiltonian of the form 



H = X] i H imftm' ~ HdC ) a tma a i'm'a (13) 
{ima} 

+ 2 y ' U mm ini m(T ni m i — a 

imm'a 

+ 2 (Umm> ~ J m m') n ima n im'a> 

im^m'a 



where H LDA is the effective Kohn-Sham-Hamiltonian derived from a self- 
consistent DFT-LDA calculation. This one-particle Hamiltonian is then cor- 
rected by Hubbard terms for direct and exchange interactions for the "corre- 
lated" orbitals, e.g. d or / orbitals. In order to avoid double counting of the 
Coulomb interactions for these orbitals, a correction term H dc is subtracted 
from the original LDA Hamiltonian. The resulting Hamiltonian (13) is then 
treated within dynamical mean field theory by assuming that the many-body 
self-energy associated with the Hubbard interaction terms can be calculated 
from a multi-band impurity model. 

This general scheme can be simplified in specific cases, e.g. in systems with 
a separation of the correlated bands from the "uncorrelated" ones, an effective 
model of the correlated bands can be constructed; symmetries of the crystal 
structure can be used to reduce the number of components of the self-energy 
etc. 

In this way, the application of DMFT to real solids crucially relies on an 
appropriate definition of the local screened Coulomb interaction U (and Hund's 
rule coupling J). DMFT then assumes that local quantities such as for example 
the local Green's function or self-energy of the solid can be calculated from a 
local impurity model, that is one can find a dynamical mean field Qq such that 
the Green's function calculated from the effective action 

S = dTj2c]na(r)Gorl mla (T -T')c m > a (T' ) 

J ° ma 

1 

+ 27 dT Umm ' 

n ma 

^ mm' a 
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(r)n mV (r) (14) 

^ JO j , 

coincides with the local Green's function of the solid. This is in analogy to 
the representability conjecture of DMFT in the model context, where one as- 
sumes that e.g. the local Green's function of a Hubbard model can be rep- 
resented as the Green's function of an appropriate impurity model with the 
same U parameter. In the case of a lattice with infinite coordination number 
it is trivially seen that this conjecture is correct, since DMFT yields the exact 
solution in this case. Also, the model context is simpler from a conceptual 
point of view, since the Hubbard U is given from the outset. For a real solid 
the situation is somewhat more complicated, since in the construction of the 
impurity model long-range Coulomb interactions are mimicked by local Hub- 
bard parameters. 1 The notion of locality is also needed for the resolution of 
the model within DMFT, which approximates the full self-energy of the model 
by a local quantity. Applications of DMFT to electronic structure calculations 
(e.g. the LDA+DMFT method) are therefore always defined within a specific 
basis set using localized basis functions. Within an LMTO [21] implemen- 
tation for example locality can naturally be defined as referring to the same 
muffin tin sphere. This amounts to defining matrix elements Gxr ^/r/ {iw) of 
the full Green's function 

G(r,x>,iu) = Yl xln(r)GLn,L>n>(iu>)xL'K'(r') (15) 

LL'RR' 

and assuming that its local, that is "on-sphere" part equals the Green's func- 
tion of the local impurity model (13). Here R, R' denote the coordinates of 
the centres of the muffin tin spheres, while r, r' can take any values. The 
index L = (n,l,m) regroups all radial and angular quantum numbers. The 
dynamical mean field Qq in (13) has to be determined in such a way that the 
Green's function Gi mpur ityL,L' of the impurity model Eq.(13) coincides with 
GL~R,L'n.(iu) if the impurity model self-energy is used as an estimate for the 
true self-energy of the solid. This self-consistency condition reads 

G impurity (i^n) = ^ (i^n + V ~ H Q (k) - S(iW n )) _1 
k 

where S, Hq and G are matrices in orbital and spin space, and iu> + ^ is a 
matrix proportional to the unit matrix in that space. 

Together with (13) this defines the DMFT equations that have to be solved 
self-consistently. Note that the main approximation of DMFT is hidden in the 
self-consistency condition where the local self-energy has been promoted to 
the full lattice self-energy. 

The representability assumption can actually be extended to other quantities 
of a solid than its local Green's function and self-energy. In "extended DMFT" 
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[22-25] e.g. a two particle correlation function is calculated and can then be 
used in order to represent the local screened Coulomb interaction W of the 
solid. This is the starting point of the "GW+DMFT" scheme described in 
section 6. 

Despite the huge progress made in the understanding of the electronic struc- 
ture of correlated materials thanks to such LDA+DMFT schemes, certain con- 
ceptual problems remain open: These are related to the choice of the Hubbard 
interaction parameters and to the double counting corrections. An a priori 
choice of which orbitals are treated as correlated and which orbitals are left 
uncorrected has to be made, and the values of U and J have to be fixed. At- 
tempts of calculating these parameters from constrained LDA techniques are 
appealing in the sense that one can avoid introducing external parameters to 
the theory, but suffer from the conceptual drawback in that screening is taken 
into account in a static manner only [26]. Finally, the double counting terms 
are necessarily ill defined due to the impossibility to single out in the LDA 
treatment contributions to the interactions stemming from specific orbitals. 
These drawbacks of LDA+DMFT provide a strong motivation to attempt the 
construction of an electronic structure method for correlated materials beyond 
combinations of LDA and DMFT. 

2. The ^-functional 

As noted in [27, 28], the free energy of a solid can be viewed as a functional 
T[G, W] of the Green's function G and the screened Coulomb interaction W. 
The functional Y can trivially be split into a Hartree part Th and a many body 
correction ^, which contains all corrections beyond the Hartree approxima- 
tion : T = Th + \f. The latter is the sum of all skeleton diagrams that are 
irreducible with respect to both, one-electron propagator and interaction lines. 
^[G, W] has the following properties: 



We present in the following a different derivation than the one given in [27]. 

We start from the Hamiltonian describing interacting electrons in an external 
(crystal) potential v ext : 



5V 

6G ~ 
SW ~~ 



(16) 




(17) 



with 




(18) 



i 
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The electron-electron interaction V ee (ri — rj) will later on be assumed to be a 
Coulomb potential. With the action 

S = - J dT&(6 T -H - VHartree)* 

+ - J dr : n (r) : V (r - r') : n (r') : (19) 

the partition function of this system reads : 

Z = J T>VB&exp(-S) (20) 

Here the double dots denote normal ordering (i.e. Hartree terms have been 
included in the first term in 18). We now do a Hubbard-Stratonovich trans- 
form, decoupling Coulomb interactions beyond Hartree by a continuous field 
cf), introduce a coupling constant a for later purposes (a = 1 corresponds to 
the original problem) and add source terms coupling to the density fluctuations 
\E'^\E' and the density of the Hubbard-Stratonovich field respectively . The free 
energy of the system is now a functional of the source fields E and P : 

F [E, P]=\nJ D^D^ t D0exp(-5[E, P]) (21) 



with 

5[E, P] = -J dT^G Hartree^ + \ J dr 

- ia J dT(j)(^y - n) + J drE^ 1 "^ + - J drPcfxj) (22) 

If in analogy to the usual fermionic Green's function G = — (T\l/\l/ ' ) we define 
the propagator W = (T<p<p) of the Hubbard-Stratonovich field <p, our specific 
choice of the coupling of the sources E and P leads to 

SF 

and 

Performing Legendre transformations with respect to G and W/2 we obtain 
the free energy as a functional of both, the fermionic and bosonic propagators 
G and W 

T[G, W] = F + tr(EG) - ^tr(PW) (25) 
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We note that W can be related to the charge-charge response function %(r, r'; r— 
r') = (r r [p(r,r)-n(r)][p(r',r')-n(r')]}: 

W(r,r';ioj n ) = V(r-r') 

- J dr 1 dr 2 V(r-r 1 )x(r 1 ,r 2 ;iuJn)V(r2-r') (26) 

This property follows directly from the above functional integral representation 
and justifies the identification of W with the screened Coulomb interaction in 
the solid. Taking advantage of the coupling constant a introduced above we 
find that T is naturally split into two parts 

F a=1 [G, W] = F a=0 + *[G,W] (27) 

where the first is just the Hartree free energy 

T a=0 (G,W) = TrlnG-TrKGx 1 -G- r )G] 

- ^TrlnW +^Tr[(v- 1 -W'^W] (28) 

(with G^ 1 = iuj n + fj, + V 2 /2 — Vh denoting the Hartree Green's function and 
Vh the Hartree potential), while the second one 

*= C ' doS- (29) 
jo oa 

contains all non-trivial many-body effects. Needless to mention, this correla- 
tion functional ^ cannot be calculated exactly, and approximations have to be 
found. 

The GW approximation consists in retaining the first order term in a only, 
thus approximating the ^-functional by 

y[G, W] = -^Tr(GWG). (30) 



We then find trivially 



<5* 

S = — = -GW (31) 



<5* 

w = Ga (32) 



GW+DMFT 



Inspired by the description of screening within the GW approximation and 
the great successes of DMFT in the description of strongly correlated materials, 
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the "GW+DMFT" method [9] is constructed to retain the advantages of both, 
GW and DMFT, without the problems associated to them separately. In the 
GW+DMFT scheme, the 'I' -functional is constructed from two elements: its 
local part is supposed to be calculated from an impurity model as in DMFT, 
while its non-local part is given by the non-local part of the GW ^-functional: 

* = *g^- loc [G RR ' , W RR '} + V imp [G RR , H^ RR ] (33) 

Since the strong correlations present in materials with partially localized d- or 
f-electrons are expected to be much stronger in their local components than in 
their non-local ones the non-local physics is assumed to be well described by 
a perturbative treatment as in GW, while local physics is described within an 
impurity model in the DMFT sense. 

More explicitly, the non-local part of the GW+DMFT ^-functional is given 

by 

^ w ~ loc [G RR ' ,W RR '} = V GW [G RR ', W RR '} - ¥£ C W {G RR ' ,W RR '} (34) 

while the local part is taken to be an impurity model ^ functional. Following 
(extended) DMFT, this onsite part of the functional is generated from a local 
quantum impurity problem (defined on a single atomic site). The expression 
for its free energy functional Ti mp [Gi mp , Wi mp ] is analogous to (27) with Q 
replacing Gh and U replacing V : 

^imp[Gi m p,W^i m p] = Trill Gi m p ~ Tr\(Q ^jmp)^'™?] 

- l -Tr In W imp + \Tr\iU' 1 - W£ p )W imp \ 

~\~ ^ ' imp\Girapi^^imp\ (35) 

The impurity quantities Gi mp , Wi mp can thus be calculated from the effective 
action: 

S = J drdr' [- £ c\{T)Qll,{r - t')c l ,{t') (36) 

+ \ E : c U T Mr) : W Ll W4 (r - t>) : c{ z {t')c u {t') :" 

where the sums run over all orbital indices L. In this expression, c^ L is a cre- 
ation operator associated with orbital Zona given sphere, and the double dots 
denote normal ordering (taking care of Hartree terms). 

The construction (33) of the ^-functional is the only ad hoc assumption 
in the GW+DMFT approach. The explicit form of the GW+DMFT equations 
follows then directly from the functional relations between the free energy, the 
Green's function, the screened Coulomb interaction etc. Taking derivatives of 
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(33) as in (15) it is seen that the complete self-energy and polarization operators 
read: 

E TC (k,iw n ) LJ/ = X x G c w (k,iu n ) LL , (37) 

" E GWr( k > ^n)w + [^ p (iuJ n )} LL , 

k 

^(q>nW = P GW (^iv n ) a( ] (38) 
q 

The meaning of (37) is transparent: the off-site part of the self-energy is 
taken from the GW approximation, whereas the onsite part is calculated to all 
orders from the dynamical impurity model. This treatment thus goes beyond 
usual E-DMFT, where the lattice self-energy and polarization are just taken to 
be their impurity counterparts. The second term in (37) substracts the onsite 
component of the GW self-energy thus avoiding double counting. At self- 
consistency this term can be rewritten as: 

E WwWlu = - £ W™l LIL , (r)G L[Ll (r) (39) 

k L X L\ 

so that it precisely substracts the contribution of the GW diagram to the impu- 
rity self-energy. Similar considerations apply to the polarization operator. 

From a technical point of view, we note that while one-particle quantities 
such as the self-energy or Green's function are represented in the localized 
basis labeled by L, two-particle quantities such as P or W are expanded in a 
two-particle basis, labeled by Greek indices in (38). We will discuss this point 
more in detail below. 

We now outline the iterative loop which determines Q and U self-consistently 
(and, eventually, the full self-energy and polarization operator): 

■ The impurity problem (36) is solved, for a given choice of Gll' an d U a p\ 
the "impurity" Green's function 

G^ p ^-(T t c l {t)c + l ,{t'))s (40) 
is calculated, together with the impurity self-energy 

S fmp = ^impl SGimp = G~ X ~ G~^ p . (41) 

The two-particle correlation function 

XL^LsU = (■ cI^cl^t) :: c\ 3 (t')c La (t') :) s (42) 
must also be evaluated. 
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■ The impurity effective interaction is constructed as follows: 

Li---L<l ■yS 

where Og iLa = (^x^x^l-B®) is the overlap matrix between two-particle 
states and products of one-particle basis functions. The polarization op- 
erator of the impurity problem is then obtained as: 

P imp = -25* imp /5W imp = U- 1 - WrJ,, (44) 

where all matrix inversions are performed in the two-particle basis B a 
(see the discussion at the end of this section). 

■ From Eqs. (37) and (38) the full k-dependent Green's function G(k, iu n ) 
and effective interaction W(q, iv n ) can be constructed. The self-consistency 
condition is obtained, as in the usual DMFT context, by requiring that 
the onsite components of these quantities coincide with Gi mp and Wi mp . 

In practice, this is done by computing the onsite quantities 

G loc {iuj n ) = ^[^(k^-S^k,^)]- 1 (45) 

k 

W loc (iis n ) = ^[V- 1 - P^ivn)]- 1 (46) 
q 

and using them to update the Weiss dynamical mean field Q and the 
impurity model interaction U according to: 

S" 1 = G Toc + ^mp (47) 

U' 1 = Wf l c + P imp (48) 

This cycle (which is summarized in Fig.l) is iterated until self-consistency 
for Q and U is obtained (as well as on G, W, T> xc and P). Eventually, self- 
consistency over the local electronic density can also be implemented, (in a 
similar way as in LDA+DMFT [29, 30]) by recalculating p(r) from the Green's 
function at the end of the convergence cycle above, and constructing an up- 
dated Hartree potential. This new density is used as an input of a new GW 
calculation, and convergence over this external loop must be reached. 

We stress that in this scheme the Hubbard interaction U is no longer an 
external parameter but is determined self-consistently. The appearance of the 
two functions U and Wi oc might appear puzzling at first sight, but has a clear 
physical interpretation: Wi oc is the fully screened Coulomb interaction, while 
U is the Coulomb interaction screened by all electronic degrees of freedom 
that are not explicitly included in the effective action. So for example onsite 
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Figure 1. Schematic representation of the inner (DMFT) self-consistency cycle of the 
GW+DMFT scheme, consisting of the construction of the impurity model effective action (36) 
from the Weiss field Q and the impurity Coulomb interaction U, the solution of the impurity 
model leading to an estimate for the impurity self-energy E imp and the polarization Pi mp and 
the self-consistency condition where local quantities of the solid are calculated and then used 
to update the impurity model. Full self-consistency of the whole scheme requires an additional 
outer cycle updating the GW Hartree Potential corresponding to the obtained electronic density 
(cf text). 
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screening is included in W\ oc but not in U. From Eq. (48) it is seen that further 
screening of U by the onsite polarization precisely leads to W\ oc . 

In the following we discuss some important issues for the implementation 
of the proposed scheme, related to the choice of the basis functions for one- 
particle and two-particle quantities. In practice the self-energy is expanded in 
some basis set {</> l} localized in a site. The polarization function on the other 
hand is expanded in a set of two-particle basis functions {^l^L'} (product ba- 
sis) since the polarization corresponds to a two-particle propagator. For exam- 
ple, when using the linear muffin-tin orbital (LMTO) band-structure method, 
the product basis consists of products of LMTO's. These product functions 
are generally linearly dependent and a new set of optimized product basis 
(OPB) [12] is constructed by forming linear combinations of product func- 
tions, eliminating the linear dependencies. We denote the OPB set by B a = 
J2ll' ^>L<t>L'^LL'- T° summarize, one-particle quantities like G and S are ex- 
panded in {4>l} whereas two-particle quantities such as P and W are expanded 
in the OPB set {B a }. It is important to note that the number of {B a } is gener- 
ally smaller than the number of {4>l4>l>} so that quantities expressed in {B a } 
can be expressed in {^l^z/}' but not vice versa. E.g. matrix elements in 
products of LMTOs can be obtained from those in the {B a } basis via the 
transformation 

W?Xl 3 u - <0M 2 l^>£> =E°U^'0^ L4 (49) 

a/3 

with the overlap matrix 01 1 l 2 = ((j)L 1 <l>L 2 \B a ), but the knowledge of the 
matrix elements Wl 1 l 2 l 3 l 4 alone does not allow to go back to the W a @. 

4. Challenges and open questions 
Global self-consistency 

As has been pointed out, the above GW+DMFT scheme involves self-con- 
sistency requirements at two levels: for a given density, that is Hartree poten- 
tial, the dynamical mean field loop detailed in section 1.1.0 must be iterated 
until Q and U are self-consistently determined. This results in a solution for 
G, W, T> xc and P corresponding to this given Hartree potential. Then, a new 
estimate for the density must be calculated from G, leading to a new estimate 
for the Hartree potential, which is reinserted into (45). This external loop is 
iterated until self-consistency over the local electronic density is reached. 

The external loop is analogous to the one performed when GW calculations 
are done self-consistently. From calculations on the homogeneous electron 
gas [31] it is known, that self-consistency at this level worsens the descrip- 
tion of spectra (and in particular washes out satellite structures) when vertex 
corrections are neglected. Since the combined GW+DMFT scheme, however, 
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includes the local vertex to all orders and only neglects non-local vertex cor- 
rections we expect the situation to be more favorable in this case. Test calcu- 
lations to validate this hypothesis are an important subject for future studies. 
First steps in this direction have been undertaken in [32]. 

The notion of locality : the choice of the basis set 

The construction of the GW+DMFT functional ^gw crucially relies on the 
notion of local and nonlocal contributions. In a solid these notions can only 
be defined by introducing a basis set of localized functions centered on the 
atomic positions. Local components are then defined to be those functions 
the arguments of which refer to the same atomic lattice site. We stress that 
in this way the concept of locality is not a pointwise (5-like) one. It merely 
means that local quantities have an expansion (15) within their atomic sphere 
the form of which is determined by the basis functions used. This feature is 
shared between LDA+DMFT and the GW+DMFT scheme, and approaches 
that could directly work in the continuum are so far not in sight. If however, 
the basis set dependence induced by this concept is of practical importance 
remains to be tested. 

Within LDA+DMFT the choice of the basis functions enters at two stages : 
First, the definition of the onsite U and the quality of the approximation consist- 
ing in the neglect of offsite interactions depends on the degree of localization 
of the basis functions. Second, the DMFT approximation promoting a local 
quantity to the full self-energy of the solid is the better justified the more lo- 
calized the chosen basis functions are. However, in the spirit of obtaining an 
accurate low-energy description one might - depending on the material under 
consideration - in some cases be led to work in a Wannier function basis incor- 
porating weak hybridisation effects of more extended states with the localized 
ones. This then leads to slightly more extended basis functions, and one has to 
compromise between maximally localized orbitals and an efficient description 
of low-energy bands. 

In GW+DMFT, some non-local corrections to the self-energy are included, 
so that the DMFT approximation should be less severe. More importantly, this 
scheme could (via the U self-consistency requirement) automatically adapt to 
more or less localized basis functions by choosing itself the appropriate U. 
Therefore, the basis set dependence is likely to be much weaker in GW+DMFT 
than in LDA+DMFT. 

Separation of correlated and uncorrelated orbitals 

As mentioned in section 1 the construction of the LDA+DMFT Hamiltonian 
requires an a priori choice of which orbitals are treated as correlated or uncor- 
related orbitals. Since in GW+DMFT the Hubbard interactions are determined 
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self-consistently it might be perceivable not to perform such a separation at 
the outset, but to rely on the self-consistency cycle to find small values for 
the interaction between itinerant (e.g. s or p) orbitals. Even if this issue would 
probably not play a role for practical calculations, it would be satisfactory from 
a conceptual point of view to be able to treat all orbitals on an equal footing. It 
is therefore also an important question for future studies. 

Dynamical impurity models 

Dynamical impurity models are hard to solve, since standard techniques 
used for the solution of static impurity models within usual DMFT, such as the 
Hirsch-Fye QMC algorithm or approximate techniques such as the iterative 
perturbation theory are not applicable. First attempts have been made in [10] 
and [40] using different auxiliary field QMC schemes and in [39] using an 
approximate "slave-rotor" scheme. These techniques, however, have so far 
only been applied in the context of model systems, and their implementation 
in a multi-band realistic calculation is at present still a challenging project. In 
the following section we therefore present a simplified static implementation 
of the GW+DMFT scheme. 

5. Static implementation 

Here, we demonstrate the feasibility and potential of the approach within 
a simplified implementation, which we apply to the electronic structure of 
Nickel. The main simplifications made are: (i) The DMFT local treatment 
is applied only to the d-orbitals, (ii) the GW calculation is done only once, in 
the form [12]: X^V = Glda ■ W[Glda], fr° m which the non-local part of 
the self-energy is obtained, (iii) we replace the dynamical impurity problem 
by its static limit, solving the impurity model (36) for a frequency-independent 
U = U(oj = 0). Instead of the Hartree Hamiltonian we start from a one- 
electron Hamiltonian in the form: H LDA - V£™ local - ±Tr£t mp (0). The 
non-local part of this Hamiltonian coincides with that of the Hartree Hamilto- 
nian while its local part is derived from LDA, with a double-counting correc- 
tion of the form proposed in [33] in the DMFT context. With this choice the 
self-consistency condition (45) reads: 

Gt oc (iu>n) = Y.\ G H~ l ^ iuJ n)-(^Gw)non-loc (50) 
k 

— (Eimp,a — -Tr CT Sj mPiO -(0) + V^° C ) ] 

We have performed finite temperature GW and LDA+DMFT calculations (within 
the LMTO-ASA[21] with 29 k-points in the irreducible Brillouin zone) for 
ferromagnetic nickel (lattice constant 6.654 a.u.), using 4s4p3d4f states, at the 
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Figure 2. Partial density of states of d-orbitals of nickel (solid [dashed] lines give the majority 
[minority] spin contribution) as obtained from the combination of GW and DMFT (see text). 
For comparison with LDA and LDA+DMFT results see [33], for experimental spectra see [34]. 

Matsubara frequencies iui n corresponding to T = 630iv~, just below the Curie 
temperature. The resulting self-energies are inserted into Eq. (50), which is 
then used to calculate a new Weiss field according to (47). The Green's func- 
tion Gf oc (T) is recalculated from the impurity effective action by QMC and 
analytically continued using the Maximum Entropy algorithm. The result- 
ing spectral function is plotted in Fig.(2). Comparison with the LDA+DMFT 
results in [33] shows that the good description of the satellite structure, ex- 
change splitting and band narrowing is indeed retained within the (simplified) 
GW+DMFT scheme. 

We have also calculated the quasiparticle band structure, from the poles of 
(50), after linearization of S(k, iuj n ) around the Fermi level 2 . Fig. (3) shows 
a comparison of GW+DMFT with the LDA and experimental band structure. 
It is seen that GW+DMFT correctly yields the bandwidth reduction compared 
to the (too large) LDA value and renormalizes the bands in a (k-dependent) 
manner. 

We now discuss further the simplifications made in our implementation. Be- 
cause of the static approximation (iii), we could not implement self-consistency 
on Wioc (Eq. (46)). We chose the value of U(u = 0) (~ 3.2eV) by calculating 
the correlation function \ and ensuring that Eq. (43) is fulfilled at lo = 0, 
given the GW value for Wi oc (lo = 0) (~ 2.2eV for Nickel [38]). This pro- 
cedure emphasizes the low-frequency, screened value, of the effective interac- 
tion. Obviously, the resulting impurity self-energy Sj mp is then much smaller 
than the local component of the GW self-energy (or than V^° c ), especially at 
high frequencies. It is thus essential to choose the second term in (37) to be 
the onsite component of the GW self-energy rather than the r.h.s of Eq. (39). 



P(E) 




20 




Figure 3. Band structure of Ni (minority and majority spin) as obtained from the lineariza- 
tion procedure of the GW+DMFT self-energy described in the text (dots) in comparison to the 
LDA band structure (dashed lines) and the experimental data of [35] (triangles down) and [34] 
(triangles up). 



For the same reason, we included V^° c in Eq.(50) (or, said differently, we im- 
plemented a mixed scheme which starts from the LDA Hamiltonian for the 
local part, and thus still involves a double-counting correction). We expect that 
these limitations can be overcome in a self-consistent implementation with a 
frequency-dependent U(u) (hence fulfilling Eq. (39)). In practice, it might be 
sufficient to replace the local part of the GW self-energy by Sj mp for correlated 
orbitals only. Alternatively, a downfolding procedure could be used. 

6. Perspectives 

In conclusion, we have reviewed a recent proposal of an ab initio dynamical 
mean field approach for calculating the electronic structure of strongly corre- 
lated materials, which combines GW and DMFT. The scheme aims at avoiding 
the conceptual problems inherent to "LDA+DMFT" methods, such as double 
counting corrections and the use of Hubbard parameters assigned to correlated 
orbitals. A full practical implementation of the GW+DMFT scheme is a major 
goal for future research, which requires further work on impurity models with 
frequency-dependent interaction parameters [41, 42, 10] as well as studies of 
various possible self-consistency schemes. 
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Notes 

1 . For a discussion of the appropriateness of local Hubbard parameters see [26] . 

2. Note however that this linearization is no longer meaningful at energies far away from the Fermi 
level. We therefore use the unrenormalized value for the quasi-particle residue for the s-band (Z s = 1). 
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